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Point source searches with the IceCube neutrino telescope have been restricted to one hemisphere, 
due to the exclusive selection of upward going events as a way of rejecting the atmospheric muon 
background. We show that the region above the horizon can be included by suppressing the back- 
ground through energy-sensitive cuts. This approach improves the sensitivity above PeV energies, 
previously not accessible for declinations of more than a few degrees below the horizon due to the 
absorption of neutrinos in Earth. We present results based on data collected with 22 strings of 
IceCube, extending its field of view and energy reach for point source searches. No significant ex- 
cess above the atmospheric background is observed in a sky scan and in tests of source candidates. 
Upper limits are reported, which for the first time cover point sources in the southern sky up to 
EeV energies. 



The spectrum of charged cosmic rays is known to 
extend to energies larger than 10 20 eV, while the sources 
and acceleration processes leading to their production 
remain unclear. Active Galactic Nuclei (AGN) are 
possible source candidates and results from the Pierre 
Auger Observatory [l[ highlight, in particular, Centau- 
rus A as an AGN near a cosmic ray hot spot in the 
southern hemisphere. Neutrinos up to EeV energies 
could be produced in hadronic acceleration processes 
at this site, predicted for example in 2]. With the 
possibility to backtrace these unique messengers, which 
propagate through space practically without deflection 
or absorption, neutrino telescopes like IceCube Q have 
the opportunity to identify such sources. 
The search for neutrino point sources relies on observing 
and reconstructing muon track directions inside the 
detector volume. Expected signals stem from charged- 
current interactions of muon neutrinos. The contribution 
from tau neutrinos is small Q and not included in the 
calculations presented here. In IceCube, Cherenkov light 
from these muon tracks is detected in digital optical 



modules (DOM) @, with 60 of these arranged vertically 
on a string and a total of 80 strings to be installed by 
2011 in the deep ice at the geographic South Pole. 
The background consists primarily of high energy muons, 
produced in extended air showers, and a much lower 
contribution of atmospheric neutrinos from the same 
interactions. A standard method to suppress the muons 
from above the horizon is a cut on the reconstructed 
track direction as in which excludes the southern 
hemisphere from IceCube neutrino point source searches. 
It also limits the observable part of the spectrum, since 
the neutrino-nucleus cross section rises with energy and 
thus leads to absorption of neutrinos going up through 
Earth with energies larger than ~ 1 PeV. To extend 
the reach of a neutrino analysis to higher energies and 
enlarge the field of view, tracks from above the horizon 
have to be included. Presented here is the first point 
source analysis which achieves this through suppression 
of lower energy background muons. This sets it apart 
from muon astronomy like in [7[ and makes this search 
sensitive to cosmic neutrinos which are expected to 
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FIG. 1: (Color online), Vp + Vp effective areas for this analysis 
with 22 strings. At positive declinations, highest energies are 
suppressed through absorption. The line for —50° to —20° is 
lowered due to strong cuts above the horizon. 
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FIG. 2: (Color online). Sensitivity to point sources with E 

spectra vs. declination, comparing this analysis with 
[| and ANTARES [U. Differences in the energy coverage at 
southern declinations are discussed in the text and Fig. [3] 



follow a harder spectrum than the atmospheric back- 
ground. A possible signal contribution from gamma-ray 
showers from above the horizon cannot be excluded [g], 
but at PeV energies this could only be expected from 
very close sources. The event selection described in the 
following was optimized for a binned study of track 
directions and requires harder cuts than the likelihood 
analysis in [f| to achieve the best sensitivity. 
Data collected with 22 strings of IceCube were analyzed 
for this work, summing to a total of 276 days livetime 
over 310 days between May 2007 and April 2008. An 
online filter provided an all-sky event sample with a rate 
of 23.6 Hz. An elaborate direction reconstruction based 
on likelihood maximization was performed offline. The 
extended probability density function used in this work 
describes the detection of first photon pulses in a DOM, 
taking into account the total number of photoelectrons 
instead of ignoring all later hits as was done in previous 
IceCube searches. Both methods are described in [9L [To| . 
The enhanced reconstruction approach leads to an 
improved angular resolution at higher energies, due 
to the importance of additional Cherenkov light from 
the increasing number of secondary particles along the 
track. The median angular resolution for the final 
event selection described below is 1.3° (1-2°) under the 
assumption of an E~ 2 (_E~ 15 ) spectrum. 
Only events with declinations > —50° were considered 
in the analysis, due to the very low signal expectations 
from the small interaction volume directly above the 
detector. The atmospheric background rises steeply 
with the angle above the horizon and is dominated by 
bundles of multiple muons. Their combined Cherenkov 
emission has a signature similar to a single very high 
energy muon track induced by a neutrino. A cut based 
on the number of photon pulses per DOM, scaled with 
the reconstructed declination and estimated energy, was 



applied to suppress this background topology. 
Further background suppression was necessary and fol- 
lowed the approach from 12j with declination dependent 
event selection criteria. Cuts on five parameters were 
simultaneously optimized for the best sensitivity to two 
template signal spectra, a standard E~ 2 and a harder 
^-1.5 -jv^g cu ^. ygjiahigg i nc l u ded two energy estimators 
which were crucial for background rejection through 
rising energy thresholds above the horizon. For upward 
going events, the efficient signal separation relied mainly 
on track quality parameters to reject misreconstructed 
muons. Final cuts were parametrized as linear functions 
of declination, smoothly connecting both hemispheres. 
The resulting sample consists of 1877 events with decli- 
nations in the range —50° to +85°. Simulations show 
that tracks coming from below are mainly induced by 
atmospheric neutrinos, while downward going events are 
dominated by atmospheric muons. The comparison of 
effective areas in Fig. [T] illustrates the fact that neutrinos 
at highest energies are absorbed prior to detection when 
coming from too far below the horizon, but are accessible 
as downward going tracks, as discussed above. 

The reconstructed event directions were kept blinded 
for the optimization of the analysis in order to avoid 
bias. In a final step, the true astronomical coordinates 
were reintroduced to determine significances for the scan 
of the sky and tests of candidate sources, described in 
the following. Details on the binned search for spatial 
clustering are given in [T3 |. The major difference is that 
in this work we extended the scan to the southern sky. 
A declination dependent optimization of the search bin 
radius showed only a weak dependence of the sensitivity 
on the exact size and led to a fixed value of 2.5°. The 
sensitivity to an E~ 2 spectrum, defined as the mean up- 
per limit at 90% confidence level constructed as in [13j, 
is displayed in Fig. [21 The steep rise in the southern 
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FIG. 3: (Color online). Left: Sky map of the directions of the 1877 events in equatorial coordinates. The color indicates the 
negative logarithm of the p-value for compatibility with background in each angular search bin before trial-corrections. Right: 
Energies for a simulated E~ 2 neutrino spectrum as a function of declination. The color indicates the fraction of events per 
angle and energy bin at final cut level. The black lines describe the region of 90% signal containment in each declination band. 



hemisphere is due to the increasing energy-sensitive cuts 
above the horizon. Signal efficiency depends strongly 
on declination and effectively defines a lower energy 
threshold rising from the TeV regime in the North 
to PeV energies in the South for an E~ 2 spectrum. 
The sensitivity can be compared to projections for 



ANTARES [111], which mainly covers GeV to TeV 
energies when looking down to the southern hemisphere. 
Since the energy range for IceCube is shifted well into 
the PeV regime above the horizon, both experiments are 
complementary in covering the full neutrino spectrum. 
No significant deviation from the background hypothesis 
was found in the sky map, shown in Fig. [3] The fluc- 
tuation with the highest significance is in the direction 
of 103.5° r.a. and 1.0° dec. (8 events observed, 1.2 
expected). By performing searches on randomized sky 
maps, we calculated the (post-trial) probability for an 
equally or more significant excess anywhere in the sky 
due to a background fluctuation to be 37%. 
We performed a separate test of pre-defined source can- 
didates, aiming at reducing the effects of trial factors. 
The selection is focused on AGN with high fluxes of 
GeV gamma-rays, motivated by hadronic models like 
L4j. By imposing thresholds for average (maximum) 



photon fluxes of 15xl0 -8 cm~ 2 s _1 (40xl0~ 8 cm^s' 1 ), 
we selected 23 bright blazars from the confirmed AGN 
in the third EGRET catalog (HI, the most complete 
list of GeV detections at the time our list was defined. 
Two additional blazars reported by Fermi in November 
2008 0, [l7| were also included based on the same 
criteria, as well as the close-by AGN M87. With the 
extension to the southern hemisphere, this is the first 
IceCube analysis that allows a search for neutrinos from 
Centaurus A and the center of our galaxy, Sgr A*, and 
both were added to the list. 

No significant excess above background was found for 



TABLE I: (color online). Results for pre-defined source can- 
didates in the southern hemisphere. /190 is the neutrino event 
upper limit of the Feldman-Cousins 90% confidence inter- 
val, $90 the resulting flux upper limit for an E~ 2 neutrino 
spectrum, i.e. d$>/dE < c&golO^GeV^cm-V^S/GeV)- 2 . 
AE is the energy range of 90% signal containment. 



Object 


r.a. [°] 


dec. [°] 


M90 


$90 


AE [GeV] 


PKS 0537-441 


84.7 


-44.1 


1.8 


281.3 


3xl0 6 


- 7xl0 8 


Centaurus A 


201.4 


-43.0 


3.7 


556.7 


2xl0 6 


- 8x10 s 


PKS 1454-354 


224.4 


-35.6 


2.0 


238.9 


1x10 s 


- 9x10* 


PKS 1622-297 


246.5 


-29.9 


4.0 


442.8 


lxlO 6 


- 9x10 s 


Sgr A* 


266.4 


-29.0 


2.3 


244.3 


1x10 s 


- 9x10 s 


PKS 1622-253 


246.4 


-25.5 


2.4 


218.3 


lxlO 6 


- 8x10 s 


PKS 1830-21 


278.4 


-21.1 


2.3 


135.5 


8xl0 5 


- 6x10 s 


PKS 1730-130 


263.3 


-13.1 


4.6 


92.4 


2xl0 5 


- 3x10* 


PKS 1510-089 


228.2 


-9.1 


3.4 


38.3 


8xl0 4 


- 2x10 s 


PKS 1406-076 


212.2 


-7.9 


4.2 


38.2 


6xl0 4 


- 2x10 s 


PKS 2023-07 


306.4 


-7.6 


2.8 


24.8 


6xl0 4 


- 2x10 s 


3C279 


194.1 


-5.8 


4.0 


27.1 


5xl0 4 


- 1x10 s 


PKS 0336-019 


54.9 


-1.8 


3.0 


14.7 


3xl0 4 


- 9xl0 7 


PKS 0420-014 


65.8 


-1.3 


3.1 


14.7 


2xl0 4 


- 8xl0 7 



any of the candidates. The upward fluctuation with 
the highest significance occurs for PKS 1622-297 with 
1 event observed and 0.3 expected. From background 
simulations, we expect a similar or more significant 
result from any of the sources in 98% of all cases. 
Flux limits, calculated according to are shown 

in Table |T] for candidates in the southern sky. They 
are considered to be valid in the given energy ranges 
containing 90% of the signal expectation. For sources in 
the northern hemisphere, the dedicated IceCube analysis 
in [6j provides better limits. The flux limit obtained for 
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Centaurus A is shown in Fig. [4] in relation to neutrino 
predictions from Cuoco & Hannestad 0] and Kachelriess 
et al. [3. 

For the blazar 3C279, we considered the possibility 
of variable emission, motivated by the reports on a 
gamma-ray flare by MAGIC Since hadronic accel- 

eration models seem to be favored for this event [20j |. 
we performed a time dependent test of neutrino fluxes 
predicted in 2l|. The most promising periods for flares 
during the livetime for the analysis were selected based 
on X-ray and optical observations from |22| , requiring 
enhanced emission in both wavebands. The chosen 
periods, given in MJDs, are 54258.6 to 54260.6, 54307 to 
54314.5, 54324.1 to 54339.9 and 54348 to 54350, adding 
up to a total of 25.4 days. To account for detector asym- 
metries and seasonal background variations affecting 
searches at timescales of tens of days or less we used a 
method of background estimation described in [23| ■ The 
strongly reduced background in the search window led 
to a rc-optimized bin radius of 4° . No events were found 
during the selected periods. We derived a limit on the 
neutrino fluence of 3.7 x 10 -3 erg cm~ 2 in the energy 
range 8.3 x 10 5 to 7.6 x 10 7 GeV, based on the model 
spectrum from [2l| shown in Fig. 21 
The searches for local excesses of events depend on 
declination dependent background expectations from 
scrambled data. Statistical errors were taken into 
account for the calculation of limits with the method 
from [241 ]. Systematic uncertainties were estimated with 
respect to the effects on the sensitivity to an E~ 2 muon 
neutrino spectrum. Variation of simulation parameters 
led to a change of —8/+ 17% for upper limit values. 
This is mainly due to remaining uncertainties about the 
light propagation in the ice as discussed in 0, [25| . To 
estimate the influence of reconstruction biases for simu- 
lated events, we compared data from below the horizon 
with a simulation of atmospheric neutrinos and derived 
an effect of — 9/+4% on the sensitivity. Together with 
theoretical uncertainties for interaction and propagation 
of muons and neutrinos (±5%), the total uncertainty 
adds up to — 13/+18% and was included in all upper 
limits and sensitivities with the method from [241 ]. 
We presented a new approach to extend neutrino point 
source searches above the horizon, which increases the 
reach of neutrino telescopes beyond PeV energies. The 
method was applied to IceCube data taken with 22 
strings. Results are consistent with the background-only 
hypothesis. The derived flux limits are the first to cover 
neutrino point sources in the southern hemisphere at 
and above PeV energies. 
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FIG. 4: (Color online). Flux limits from this work for 
2 sources: Centaurus A limit (blue dashed) with models 
CH08 @] and KOT08 QJ]; 3C279 flare search (red solid) with 
model Reimer09 [H. 
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